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Recent experiments have probed quantum dots through transport measurements in the regime 
where they are described by a two lead Anderson model. In this paper we develop a new method to 
analytically compute for the first time the corresponding transport properties. This is done by using 
the exact solvability of the Anderson Hamiltonian, together with a generalization of the Landauer- 
Biittiker approach to integrable systems. The latter requires proper identification of scattering 
states, a complex and crucial step in our approach. In the Kondo regime, our results include the 
zero-field, finite temperature linear response conductance, as well as the zero-temperature, non- 
equilibrium conductance in an applied Zeeman field. 



In recent years there have been a flurry of experimen- 
tal studies of quantum dots [Qj2) . In these experiments a 
single localized dot state interacts with connecting leads 
notwithstanding that a finite number of electrons sits on 
the dot. In a testament to advances in the miniatur- 
ization of solid state technology, the occupancy of this 
level can be controlled by a gate voltage, thereby con- 
trolling the current through this 'single-electron transis- 
tor'. For an odd number of electrons, the dot becomes, 
due to Kramers degeneracy, a nanoscale realization of 
the single-impurity Kondo system - one has in effect a 
'tunable' Kondo model The observed split peak 

structure in the non-equilibrium conductance in an ap- 
plied Zeeman field is a hallmark of the Kondo physics 

The experimental probes of these systems are measure- 
ments of transport properties. The appropriate theoret- 
ical framework to describe electronic transport through 
the dot is the two-lead Anderson model (to distinguish 
it from the one-lead model appropriate for studying di- 
lute impurities in a bulk metal). This model has been 
examined using a variety of techniques. In the T = 
linear response conductance was studied both in and out 
of a magnetic field via the Friedel sum rule (although the 
occupancy of the dot necessary to apply the rule was com- 
puted only in perturbation theory). Using the NCA ap- 
proximation and perturbation theory, studied quan- 
tum dots out-of-equilibrium. However such techniques 
have limitations: the NCA approximation does not work 
in the presence of a magnetic field while perturbation 
theory requires either the dot Coulomb repulsion or the 
dot-lead coupling to be small. Numerically, the zero- 
field, finite temperature linear response conductance has 
been computed in using the numerical renormaliza- 
tion group (NRG). However the latent integrability of 
the (one-lead) Anderson model Jh],[ii[] has never been 
exploited to compute directly transport properties. 

In this work we combine the integrability of the An- 
derson model with Landauer-Buttiker transport theory. 



To do so one must first face the challenge of identifying 
scattering states in the context of integrability. In |8|] 
the Bethe ansatz of the simpler one-lead Kondo model 
was employed to study the equilibrium impurity DOS, of 
which the non-equilibrium counterpart, in context, would 
yield the desired transport properties j^]. Neither the 
computation of the latter nor the scattering states in the 
two lead model were addressed in The identification 
of the scattering states is non-trivial. The correct scat- 
tering states will necessarily be electronic in nature (i.e. 
carry charge e and spin 1/2) and will be confined to a 
single lead. In contrast, the eigenstates of the exact so- 
lution affect a spin-charge separation (i.e. they are not 
electronic) and they are not confined to a single lead. 
One purpose of this letter is to outline how one can un- 
derstand scattering states in terms of the eigenstates of 
the Bet he ansatz. 

We are successful in this endeavour on two counts. 
Firstly, we are able to construct scattering states at the 
Fermi surface for a dot under arbitrary gate voltage and 
arbitrary magnetic field and so are able to offer for the 
first time a proof of the Friedel sum rule based upon in- 
tegrability, alternative to that from Fermi liquid theory 
§. Secondly, provided the gate voltage is adjusted so 
that the dot sits at the symmetric point, we can compute 
scattering states away from the Fermi surface. Through 
a Landauer-Buttiker approach we are then able to com- 
pute the zero-field finite temperature linear response con- 
ductance (in agreement with numerical results of 0), as 
well as the T = non-equilibrium magneto-conductance, 
exhibiting the observed peak splitting. 

Model: To then begin, the two-lead Anderson model 
Hamiltonian is given in the continuum limit by 

/oo 
dx{-ic\ ry {x)d x ci IJ (x) + ViS(x)[cj cr (x)d a 
la -°° 

+ rft Cfg- (a:)] } + e d }^ n a + Un^ni, (1) 
where n a = dJ.d CT . The c/'s are the lead electrons while 
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the d's are the dot electrons. Here J^i is a sum over 
the two leads (I = 1,2). Our formalism allows for the 
possibility that the hopping matrix element, Vi, differs 
between the leads (as is typical in any experimental re- 
alization). For simplicity, in our presented results we 
assume V\ — V2. Rather than treating the leads as half- 
lines with both left and right moving fermions, we repre- 
sent the leads as 'unfolded' with fermions that are solely 
right-moving. Fermions in cither lead that are incident 
upon the dot are considered to lie in the region, x < 0, 
while those traveling away from the dot in either lead are 
found at x > 0. We represent this in Figure 1 by drawing 
the leads as elongated arcs. 




x=0 



FIG. 1. A sketch of two leads attached to a quantum dot. 

It will be advantageous to reformulate this problem as 
a one-lead Anderson Hamiltonian. To do so, we introduce 
even/odd electrons: c e / = (V1C11LV2C2)/ ' \/V 2 + V 2 ■ Re- 
casting H. in this new basis, the odd electron, c , decou- 
ples and we are left with 

H = J2 I dx{-ict a (x)d x c ea (x) + (V 2 + Viy'Hix) 

x [cl a (x)d a + dlc ea {x)] } + e d ^ n a + Un r n L . (2) 

a 

We have thus reduced the problem to that studied using 
Bethe ansatz in a series of papers by Kawakami and Okiji 
and Filyov, Wie gmann and Tsvelick 0. 

As such we summarize briefly the results of this work. 
Applying the Bethe ansatz to a system with N particles 
and total spin 2S Z = (N — 2M) yields a set of quantiza- 
tion conditions describing a finite number of bare excita- 
tions in the system: 

AkjL+iSlkj) = TT flOj) - K+i/2 _ 
jJiS(%)-A Q -V2' 

-pr X a - g{k j ) + i/2 = _ TT X a - Xp+i 

11 X a - g{ki) - i/2 1 = 1 X a — Xp — i 1 1 > 

where S(k) = -2tan- 1 (r/(fc - e d )), g{k) = (k - e d - 
U/2) 2 /2UT, and T = {V? + V 2 2 ). 

When ed > —U/2 (the case ed < —U/2 can be handled 
through particle-hole transformations), the zero temper- 
ature ground state of the system is formed from N — 2M 
real k/s and M real A Q 's. Associated with each A Q is a 
pair of complex k's, fc^, related to the A Q via g(k±) = 

g(x(X a ) t iy{X Q )), x(X) = U/2 + e d - VUT{X + (A 2 + 
l/4) 1 /2)i/2 ; and y(A ) = y/uf(—\ + (A 2 + l/4) 1 / 2 ) 1 /2. 
Roughly speaking, the real fc's can be thought of as 
charge excitations and the A a 's as spin excitations cou- 
pled to charge excitations, i.e. bound states of electrons. 



The finite temperature ground state is considerably more 
complicated, involving an infinite hierarchy of excitations 
as categorized in the 'string-hypothesis' JTo| , pT| . 

Scattering States: Under this one-lead reformulation, 
we are still able to make contact with the scattering am- 
plitudes of electronic excitations off the quantum dot. 
Let T(e)/R(e) be the transmission/reflection amplitudes 
of electronic excitations of energy e between leads in 
the original two lead picture. On the other hand, the 
even/odd excitations will scatter off the dot with some 
pure phase, S e (e)/S (e), where in particular 5 (e) = 0. 
The two sets of amplitudes are related straightforwardly: 

e l5 ° {c) = l = R{e)-T{e). (4) 

We will exploit extensively the fact that the determi- 
nation of the phase, 5 ei gives T and R in the original 
problem. 

To determine <5 e (e), we employ an energetics argument 
of the sort used by N. Andrei in the computation of the 
T = magnetoresistance of impurities in a bulk metal 
p3[ . Imagine adding an electron to the system. Through 
periodic boundary conditions, its momentum is quan- 
tized, p = 2-nn/L. If the dot was absent, the quanti- 
zation condition would be determined solely by the con- 
ditions in the bulk of the system and we would write, 
Pbuik = 2nn/L. Upon including the dot, this bulk mo- 
mentum is shifted by a term scaling as 1/L. The quan- 
tization condition is then rewritten as 

p = 2im/L = p bulk + S e (e)/L, (5) 

where L is the system's length. The coefficient of the 
1/L term is identified with the scattering phase of the 
electron off the dot. 

In order to determine the scattering phase of an elec- 
tron (as opposed to a spin or charge excitation), we must 
must specify how to glue together a spin and a charge ex- 
citation to form the electron. The situation is analogous 
to adding a single particle excitation in the attractive 
Hubbard model. Adding a single spin f electron to the 
system demands that we add a real k (charge) excita- 
tion. But at the same time we create a hole at some 
A in the spin distribution. The number of the available 
slots in the spin distribution is determined by the num- 
ber of electrons in the system. Adding an electron to 
the system thus opens up an additional slot in the A- 
distribution. The electron scattering phase off the impu- 
rity is then the difference of the right-moving k-impurity 
momentum, p- lmp (k), and the left-moving A-hole impurity 
momentum — Pi mp (A): 

5 l = Pimp = Pimp(^) + Pimp (A). (6) 

We now turn our attention to computing these impurity 
momenta. As part of this, we will relate the impurity mo- 
menta to the impurity density of states (which in turn, 
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will allow us to prove the Friedel sum rule). We at the 
same time will compute the energy of the k and A ex- 
citations in order to parameterize scattering in terms of 
energy. 

In the continuum limit, these excitations are described 
by smooth densities, p(k) for the real k's and cr(A) for 
the A's. From (||), equations describing these densities 
can be derived in the standard fashion loJlOl : 



P(k) = ^ + ^ + g'(k) J dX a x {g{k) - X)a(X); 
a(A)= _^W + ^)_ f Q dX'a 2 (X' - A)er(A') 



dk ai(A - g(k))p(k), 



(7) 
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where L is the system size and A(k) — dk5(k)/2ir, 
A(A) = — d\Ke5(x(X) + iy(X))/n, and a n (x) = 
2nj (7r(n 2 +4x 2 )). B and Q mark out the 'Fermi-surfaces' 
of the k and A distributions. A key observation to make 
of these equations is that one can divide the densities into 
bulk and impurity pieces via p(k) — > Pbu\k(k) + Pi mp (k) / L 
and similarly for a(X). The impurity densities of states 
contain all the information needed about degrees of free- 
dom living on the quantum dot. For example the to- 
tal numbers of spin | and J, electrons living on the dot 
are n dT = Jq dXa imp (X) + J B Qo dkp imp (k) and n dl = 

Iq dXai mp (X). 

The energies and momenta of these excitations can be 
derived through well known techniques |l2| . The energies 
are given by 



e(k) = k 



H 



dXe(X) ai (X - g(k)); 



Q 



e{X) = 2x{X)- / e(A> 2 (A'-A) 

B 

g'(k)e(k)ax(g(k) — A). 



(8) 



The momenta are akin to the densities in that they divide 
into bulk and impurity pieces. The impurity momenta, 
the momenta relevant to scattering, are given by |l2[ 



Wmp(fc) = 6(k) + / dX a imp (X)6 1 (g(k) - A); 

Jq 

f Q 

p imp {X)=2Re6(x(X)+iy(X)) + / dX'<j imp (X')9 2 (X-X r ) 

Jq 

+ dk /9i 1 „ P (fc)6' 1 (A- 5 (fc)), 
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where 9 n (x) = 2 tan _1 (a;/n) — 2ir. 

We now observe a relationship between the impurity 
DOS and the impurity momenta key to the actual com- 
putation of the scattering phases: 



dkPim P (k) = 2irp imp (k); 
d\Pi mp (X) = -27rcr imp (A). 



(9) 



The scattering phase for a spin | excitation created from 
a k-particle and A-hole is then given to be 

pk />oo 

5\=2tt dk Pimp {k) + 2ir dX'a imp (X'). (10) 



If k and A are chosen to be at the Fermi surface, i.e. k = 
B and A = Q, we prove the Friedel sum rule, i.e. 5\ = 
2i:nd\- Through a particle-hole transformation, we can 
similarly characterize the scattering of spin j excitations 
and so also prove the Friedel sum rule in this case. 

To discuss scattering away from the Fermi surface, re- 
consider (|Io[). Suppose we want to compute the scatter- 
ing of a spin | electron with energy, e e ;. If we choose 
the pair, (fc, A), such that e(A) + e(fc) = e e i, we will have, 
via ( |To| ) , the scattering of an excitation with energy, e e i . 
However we encounter a problem in that there is not a 
unique pair, (/c,A), for a given energy. While we cannot 
solve the problem of choosing the correct pair (or com- 
bination of pairs) in general, we can make progress when 
we are in the Kondo regime of the Anderson model (i.e. 
U + 2ed ~ 0). In this regime we expect the scattering 
phase to vary on the scale of the Kondo temperature, 
Tfc. The electron scattering phase is determined by yOi mp 
and <7i mp , the two impurity densities. Of the two, only 
Pimp varies on scales on the order of 71- (eri mp is con- 
trolled by the much larger scale, y/UT, governing charge 
fluctuations). Thus in computing electronic scattering 
phases away from the Fermi surface at zero temperature, 
it is natural to keep A = Q, its Fermi surface value, and 
vary k. Specifically, to describe an electron with energy, 
e e i, we choose (k,X) such that A = Q and k such that 
e(k) = e e i. With this ansatz, we then have restricted the 
two dimensional phase space, (A, fc), of potential excita- 
tions carrying the quantum numbers of an electron to an 
one dimensional subspace. We have further arguments 
that suggest this ansatz is exact at the symmetric point 
of the Anderson model |[2| . 

Linear Response (T > 0): We now try out this ansatz 
by examining a quantity that requires us to understand 
scattering at finite energy: the linear response conduc- 
tance as a function of temperature. We compare it to 
Costi et al.'s NRG results and find excellent agree- 
ment. This is important as it indicates we have an es- 
sentially correct description of the low energy scattering 
states. 

Computing the linear response conductance at finite T 
is a complicated matter. We now have to compute the 
scattering phases of the glued charge and spin excitations 
in the presence of a 'thermalized ground state'; that is we 
must compute dressed scattering matrices. This ground 
state is no longer composed of merely real k states and 
A states of bound spin and charge as it was at T = 0. 
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Rather all the possible solutions of the Bethe ansatz 
equations of the model make an appearance. 
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FIG. 2. A plot of the scaling curve for the conductance at 
the symmetric point as a function of T/Tk- 71- is as defined in 
Costi et al.'s work and as such there are no free parameters. 

However it is possible to derive equations governing 
the impurity densities p lmp and u luip in the presence of 
the complicated ground state and to solve numerically 
the corresponding set of coupled integral equations [Q . 
With these in hand, the linear response conductance is 
given by 

/oo 
de el (-d eel f(e el ))\T(e el )\ 2 , (11) 
-oo 

where / is the Fermi distribution and |T(e e ;)| 2 = 
sin 2 (i<5 e z(e e 2, T)) is the dressed scattering amplitude. A 
plot of the result and a comparison with the NRG re- 
sults of § are given in Figure 2. For temperatures up to 
Tk, the regime where in principle an NRG computation 
should be of greatest accuracy, we find excellent agree- 
ment. 

Because of the Fermi liquid nature of the problem when 
T<T(., we know the functional form of the conductance 
is 

G(T/T k ) = 1 - cT 2 /T k 2 + ■■■. (12) 

Costi et al. , based upon results borrowed from @,|l5| , 
computed c — 7r 4 /16 — 6.088. We find from a fit of our 
curve, c = 6.05 ± .1. We have arrived at this value by 
fitting the plot in the region T/T k < .1. The error is 
systematic in nature, arising from the arbitrary nature 
of deciding the region over which to fit. 

We also compare our results in Fig. 2. withp7). It 
would appear that the logarithmic dependence ^[char- 
acteristic of weak coupling and arising from a one-loop 
RG, should only be expected to become qualitatively de- 
scriptive for values of T/T k in excess of about 20. 

We have recently computed the finite temperature lin- 
ear response conductance at the symmetric point in finite 
applied field. We again find reasonable agreement with 

0- 



N on- Equilibrium (T = 0, H ^ 0): As we have suc- 
cessfully reproduced numerics on the finite temperature 
linear response conductance, we have some confidence 
that we understand scattering away from the Fermi sur- 
face. As such we can explore non-equilibrium conduc- 
tances which depend upon the same information. Al- 
though there are issues relating on how to view the non- 
equilibrium system in the one-lead picture, these can 
be successfully handled fj^ |. We essentially adopt the 
Landauer-Biittiker approach taken in Jl6[ : we employ 
the in-equilibrium scattering matrices; the sole role the 
bias plays is to set the particle distribution in each lead. 
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FIG. 3. A plot of the differential conductance in a magnetic 
field at the symmetric point with U = 871T. 

The current between the two leads, lead one at zero 
bias and lead two at bias V < is given by 

J(V) = lj\e{\T^ 2 {e)\ 2 + |T^ 2 ( e )| 2 ), (13) 

In Figure 3 we plot the differential conductance, G — 
—dyJ, in the presence of a magnetic field. Our results 
are in rough accordance with § - we find that for fields, 
H > Tk, the zero-bias, zero field differential conduc- 
tance peak divides in two, one peak for each spin species. 
Roughly speaking, the origin of the split in the differen- 
tial conductance arises from a similar bifurcation in the 
impurity density of states. The spectral weight of the 
Kondo resonance present at zero energy when H = 
divides into two resonances near V ~ ±iJ , again one as- 
sociated with each spin species. However unlike ^ we 
find that the peak does not occur exactly at \V\ = H - in 
fact, the approach to this value is logarithmic in H [12[ . 
This discrepancy is unsurprising given the perturbative 
nature of the computation in [^) . 

In the large field limit, H ^> Tfc, our analytical con- 
trol allows us to employ a Wiener-Hopf analysis [jl2| to 
analyze the properties of the peak in the differential con- 
ductance. With a = log(i7v / 7re/2 3 / 2 Tfe)/7r, a peak max- 
imum occurs at a bias, eV maK = —H(l — cot _1 (a)/(27r)). 
The corresponding width of the peak is given by eAV = 
g/(27r)(cot~ 1 (a-l/2)-cot' 1 (q + l /2)), wh ile the height 
of peak is, G max = e 2 //i(3/2 - a/V4a 2 + 1). 
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To conclude, we have shown that the integr ability of 
the Anderson model can be exploited to compute trans- 
port quantities for quantum dots. The computation is 
more difficult than in the quantum Hall tunneling prob- 
lem [ [l6[ , and requires an ansatz to describe scattering 
away from the Fermi surface. The accuracy obtained in 
reproducing known numerical results is remarkable, 
and suggests that the method could be a successful way 
of handling non-equilibrium transport in general inte- 
grable systems. In future work, we will further exploit 
this technique to compute the DC shot noise at zero tem- 
perature in quantum dot systems Q . 
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